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一 种 用 于 空间 VLBI 的 射电 源 条 纹 搜 索 算法 


ESR, 童 Fess HAR, ak 娟 5245 
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京 210023; 5. 上 海 市 导航 定位 重点 实验 室 ， 上 海 200030) 

摘要 : 我 国 探 月 工程 嫦娥 七 号 将 构建 由 中 继 星 上 4. lmD 径 抛物 面 天 线 及 地 基 射 电 望远镜 
组 成 的 首 个 月 地 空间 VLBI 试 验 系统 。 不 同 于 地 基 VLBI 望 远 镜 , 受 中 继 星 轨道 扰动 和 设备 时 延 
影响 ， 射 电源 预报 时 延 模型 和 实际 时 延 有 时 可 能 会 有 较 大 差异 ,导致 无 法 引导 VLBI 相 关 处 理 
机 正常 工作 ,需要 通过 射电 源 条 纹 搜 索 方 法 寻找 实际 时 延 模型 才能 解决 这 一 问题 . 由 于 4. lm 
空间 天 线 口 径 较 小 ， 接 收 信号 能 力 弱 。 为 提高 信 骂 比 ， 需要 长 时 间 积 分 ， 而 线性 时 延 模型 无 
法 保证 长 时 间 的 时 延 精度 。 针 对 这 一 问题 , 本 文 提出 了 一 种 用 于 空间 VLBI 的 射电 源 条 纹 搜 索 
算法 ， 首 先 通过 建立 时 延 和 时 延 率 的 二 维 搜索 网 格 ,使 用 每 个 网 格 点 对 应 的 先 验 时 延 和 时 延 
率 构 建 线 性 时 延 模 型 ; 再 使 用 该 线性 时 延 模型 对 原始 数据 分 时 段 相 关 处 理 ; 然后 采用 二 维 傅 
里 叶 变 换 搜 索 残 余 时 延 和 时 延 率 ,确定 最 大 的 相关 幅度 值 对 应 的 网 格 点 ; 最 后 使 用 该 网 格 点 
的 先 验 时 延 和 时 延 率 以 及 残余 时 延 和 时 延 牵 重 构 二 次 项 时 延 模 型 。 通 过 对 RadioAstron 空 间 
VLBI 的 地 空 基线 实际 数据 处 理 ， 验 证 了 该 射电 源 条 纹 搜 索 算 法 的 正确 性 。 
关键 字 空间 VLB1; 条 纹 搜索 ， 时 延 模型 ， 相 关 处 理 ， 射 电源 
中 图 分 类 号 : TN216 ”文献 标识 码 : A 文章 编号 : 1007-2276- (2004) 4-0338-05 
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其 长 基线 干涉 测量 (Very Long Baseline Interferometry, VLBI) 是 一 种 具有 高 分 辨 
率 、 高 精度 特点 的 射电 干涉 技术 "。 它 采用 原子 钟 控制 的 高 稳定 独立 本 振 系 统 和 数据 采集 装 
A, 通过 多 个 望远镜 同时 接收 目标 源 信和 号， 记录 数据 并 发 送 到 异地 的 VLBI 数 据 处 理 中 心 , 经 
处 理 得 到 观测 结果 。 在 VLBI 观 测 中 ,两 台 射 电 望 远 镜 间 的 距离 称 为 基线 ， 它 决定 VLBI 的 最 大 
分 辩 率 ”， 若 有 空间 望远镜 与 地 面 望 远 镜 构成 空地 YLBI 系 统 ， 基 线 长 度 超过 地 球 直径 ， 可 实 
现 更 高 的 分 辩 率 。 我 国 探 月 工程 四 期 任务 将 利用 嫦娥 七 号 中 继 星 上 的 4. lm 口径 抛物 面 天 线 ， 
进行 月 地 空间 VLBI 试 验 。 不 同 于 常规 地 基 YVLBI™”"， 受 中 继 星 轨道 扰动 和 设备 时 延 影响 ， 射 电 
源 预报 时 延 模型 可 能 会 和 实际 时 延 有 较 大 差异 , 无 法 引导 相关 处 理 机 正常 工作 。 通过 射电 源 
条 纹 搜索 , 可 确定 实际 的 相关 处 理 机 所 需 的 高 精度 时 延 模型 ， 用 于 空间 VLBI 空 地 基线 相关 处 
理 。 


目前 , 在 探 月 工程 中 ,进行 探 测 器 信号 相关 处 理 时 , 采用 DOR (Differential One-way 
Ranging) 信 号 的 点 频 信号 特征 来 设计 条 纹 搜索 算法 ”。 而 射电 源 信 号 是 白 噪声 信号 ， 无 法 
利用 现 有 探测 器 信号 搜索 算法 。 目 前 针对 射电 源 信 号 的 条 纹 搜索 有 eFFT 条 纹 搜 索 算法 、CAF-W 
条 纹 搜索 算法 ”和 多 重 网 格 条 纹 搜索 算法 “。eFFT 条 纹 搜索 算法 和 CAF-W 条 纹 搜索 算法 可 实 
现 残 余 时 延 和 时 延 率 的 计算 , 但 主要 用 于 强 源 的 时 延 模 型 搜索 。 多重 网 格 条 纹 搜 索 算法 的 基 
本 思想 是 建立 多 层次 、 不 同 尺 寸 的 网 格 ， 通 过 这 些 网 格 ， 在 不 同 范围 内 搜索 需要 的 目标 值 。 
多 重 网 格 条 纹 搜索 算法 使 用 线性 时 延 模型 进行 相关 处 理 ， 在 时 延 - 时 延 率 平面 中 搜索 条 纹 ， 
适用 于 信和 号 较 强 , 短 时 间 积 分 就 能 得 到 条 纹 的 应 用 场景 由 于 空间 小 口径 天 线 接收 信号 较 弱 ， 
想 要 获得 干涉 条 纹 需要 进行 较 长 积分 时 间 。 使 用 通常 的 线性 时 延 模型 无 法 保证 较 长 时 间 范 围 
内 的 相关 处 理 机 模型 时 延 精度 , 需 设法 构造 精度 更 高 的 二 次 项 时 延 模型 。 如果 在 常规 多 重 网 
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格 中 增加 时 延 二 次 项 ， 会 使 搜索 空间 从 二 维 扩展 到 三 维 ， 计 算 复 杂 度 增加 一 个 量 级 。 

本 文 基于 多 重 网 格 条 纹 搜 索 算法 , 提出 了 一 种 用 于 空间 YLBI 的 射电 源 条 纹 搜 索 算 法 : 首 
先 建立 关于 时 延 和 时 延 率 的 二 维 搜索 网 格 , 使 用 每 个 网 格 点 对 应 的 先 验 时 延 和 时 延 率 信息 构 
建 线性 时 延 模型 , 对 原始 数据 分 时 段 相关 人 处理; 然后 由 二 维 傅 里 叶 变 换 求 解 残 余 时 延 和 时 延 
率 后 进行 补偿 ,得 到 每 个 时 段 的 时 延 和 时 延 率 ,计算 每 个 网 格 点 对 应 的 相关 幅度 值 ， 最 后 在 
最 大 相关 幅度 值 对 应 的 网 格 点 上 重 构 二 次 项 时 延 模型 。 采 用 该 算法 ， 对 俄罗斯 RadioAstron 
空间 VLBI 的 实际 观测 数据 进行 相关 处 理 时 取得 了 较 好 的 效果 。 


1 射电 源 条 纹 搜索 算法 


目前 中 国 VLBI 网 软件 相关 机 使 用 五 次 项 时 延 模型 进行 相关 处 理 。 在 条 纹 搜索 的 过 程 中 ， 
为 减少 计算 复杂 度 , 通常 使 用 由 时 延 常 数 项 和 时 延 一 次 项 系数 构成 的 线性 时 延 模型 。 但 随 着 
时 间 范 围 的 扩大 ,线性 时 延 模型 的 误差 与 五 次 多 项 式 模型 的 误差 也 会 增加 。 综合 考虑 到 条 纹 
搜索 计算 量 与 模型 精度 因素 , 我 们 将 设法 构建 由 常数 项 、 一 次 项 系数 和 二 次 项 系数 构成 的 二 
次 项 时 延 模型 。 
1. 1 相关 处 理 算 法 

假设 两 全 站 接收 到 的 目标 源 射 频 信号 分 别 为 s1(t)，sz(t)， 同 一 波 前 信号 到 达 两 台 站 时 
BEF AT, M) 


S2(t) = Ss1(t + 7) (1) 

由 射频 信和 号 得 到 的 基 频 信和 号 
x (t) = s1(t)e ht (2) 
Xa (t) = sa (t)e 2" = s (t + Te 20 (3) 


HP, PAWEWRTEME. 
相关 处 理 以 台 站 1 为 参考 ， 使 用 预报 模型 ze 对 台 站 2 进行 时 延 补偿 ， 即 


Xa (t — Te) = S3 (t — Te)e or = Si 人 十 工 一 To)e-72m1otrro) (4) 
经 过 条 纹 旋转 后 : 
Xs(t — To)e hore = s(t— Tt)e iot = s(t+t— Tt)e ht (5) 


对 台 站 1 和 人 台 站 2 的 信号 进行 傅 里 叶 变 换 ; 


AOSE 
xa(t — Tee rhre 3 S, (f + fo)e PMU tote = S (f + foje TOt- (6) 
Sif + 万) 为 台 站 1 信号 经 过 频率 转换 后 的 功率 谱 , KIARA Bah EK AR N : 
Pa(P) = |Si(f + f | e+ (7) 


由 式 (7) 可 知 ， 当 残余 时 延 re 一 ft 达 到 最 小 值 时 ， 两 人 台 站 的 互相 关 功 率 谱 达 到 最 大 峰 
值 。 
1. 2 残余 时 延 搜索 算法 
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理 后 的 互相 关 功 率 谱 。 
KAT (8) 离散 化 ， 即 


F(a,,b,) =| 
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1 Residual dolay 


图 1 2-D 全 里 叶 变 换 残 余 时 延 搜索 图 

Fig. 1 Search function 

根据 二 维 傅 里 叶 变换 的 性 质 可 知 ， 对 相关 处 理 后 得 到 的 互相 关 谱 进行 二 维 传 里 时 变换 
后 ， 可 得 到 一 个 关于 残余 时 延 和 时 延 率 的 二 维 互相 关 函 数 ”“， 如 图 1 所 示 ， 其 公式 如 下 : 


POT At) = [JP {J Prof ge ntd eeir dy | 


HF, At ARRAY RE, ATAGRRIN KER, BAW TE, 7 为 一 组 二 维 傅 里 时 变换 的 数据 时 
长 ， 态 为 天 空 频率 ，F(Ar A1i) 为 不 同 残 余 时 延 和 时 延 率 下 的 相关 幅度 值 ，P12(f,t) 为 相关 处 


2BT es = 
ADE PE 


P > (n, m)e ~’? brmT/M) e-i4narnB/N | 


RP, ap Ally DAA ABEE IRR SEA IY SE 8, WE EI AR A IA ER 
傅 里 叶 变 换 点 数 。F(ay,by) 为 经 过 二 维 傅 


1.3 条 纹 搜索 与 二 次 项 时 延 模型 台 


EB 叶 变 换 后 得 到 的 相关 幅度 值 ”。 


E 构 流 各 


根据 以 上 原理 ， 对 于 空间 VLBI 射 电源 条 纹 搜索 与 二 次 项 时 延 模型 重 构 流 程 如 图 2 所 示 。 
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Constructing 2D search grid 
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图 2 射电 源 条 纹 搜索 搜索 与 二 次 项 时 延 模型 重 构 流 程 图 
Fig.2 Flow chart of radio source fringe search and quadratic term delay model reconstruction 
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4B’ 4B 2foT’ 2foT 
实测 数据 的 处 理 经验 ， 设 置 时 延 步 长 为 
dstep = 08 
drstep = 


条 纹 搜索 的 步 长 与 时 延 模型 精度 和 计算 量 密切 相关 ， 如 果 搜 索 步 长 设置 过 大 ， 可 能 得 
到 的 模型 精度 太 低 无 法 用 于 后 续 数据 处 理 ; 如 果 步 长 设置 太 小 ,由 于 目前 的 算法 采用 全 网 格 
搜索 策略 ， 依 次 计算 每 个 网 格 点 数据 会 极 大 影响 整体 的 计算 效率 。 注意 到 步 长 设置 较 小 时 可 
能 会 有 多 组 网 格 点 搜索 到 满足 精度 要 求 的 模型 后 续 会 研究 基于 二 维 随机 游 走 的 网 格 搜索 策 
略 ， 当 找到 合适 的 时 延 模型 后 就 停止 计算 ， 这 样 可 以 有 效 的 提升 效率 。 

确定 步 长 后 ， 根 据 搜索 范围 和 先 验 时 延 模型 划分 网 格 点 。 每 个 网 格 点 对 应 不 同时 延 和 
时 延 紊 ， 代 表 经 过 不 同时 延 补偿 的 数据 处 理 。 如 图 3 所 示 ， 以 先 验 时 延 模 型 中 的 时 延 和 时 延 
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率 作为 中 心 点 ， 等 步 长 计算 每 个 网 格 点 对 应 的 时 延 和 时 延 率 值 。 
由 于 弱 信 号 使 用 线性 时 延 模型 进行 长 时 间 积分 无 法 保证 时 延 精度 ， 获 得 干涉 条 纹 ， 所 
以 将 原始 数据 划 为 7 组 较 短 的 数据 段 进行 处 理 ， 每 一 组 为 一 个 二 维 傅 里 叶 变 换 组 ， 时 间 长 度 
为 7。 在 每 个 网 格 点 上 构造 线性 时 延 对 原始 数据 补偿 后 进行 相关 处 理 ， 相 关 处 理 的 积分 时 间 
为 T/M， 一 个 二 维 傅 里 叶 变 换 组 的 数据 大 小 为 M x NW。 得 到 两 台 站 的 互相 关 功 率 谱 后 ， 利 用 
TK (9) 对 每 一 组 的 相关 谱 进 行 二 维 傅 里 叶 变 换 ， 得 到 每 个 网 格 点 上 的 I 组 残余 时 延 、 时 延 率 
和 相关 幅度 值 , 利用 最 大 相关 幅度 的 平均 值 确定 最 优 解 对 应 的 网 格 点 , 在 该 网 格 点 上 进行 时 
延 和 时 延 率 的 拟 合 ， 得 到 最 优 的 二 次 项 时 延 系数 。 


> @ A and point 


图 3 二 维 网 格 示意 图 
Fig.3 Schematic diagram of two-dimensional grids 
在 拟 合 时 ， 将 每 一 组 的 残余 时 延 af AEN ESP, GRADE ADE POR NY EA A ee EN 
延 af 和 时 延 率 bf 相 加 后 ， 可 得 到 该 组 的 时 延 g; 和 时 延 率 b:， 即 
a; = af + aj (12) 
b; = bf + b? (13) 
VLBI 观 测 具 有 时 间 连 续 性 ， 时 延 模 型 也 应 该 是 连续 的 。 但 在 求解 残余 时 延 的 过 程 中 ， 
数据 分 段 导 致 每 一 组 之 间 的 时 延 具有 间断 点 , 在 拟 合 二 次 项 时 延 系数 时 , 需要 根据 连续 的 约 
束 条 件 ， 即 保持 第 一 组 的 时 延 ao 不 变 ,修正 其 余 组 的 时 延 Qi; 使 之 与 前 一 组 末 的 时 延 相等 ， 如 
式 (14) 所 示 。 


mM — 
{ ao = ao 0<i<l (id) 


al” = ai1 + biT 
其 中 ，a 了 ?为 改正 后 的 时 延 ，1 为 二 维 傅 里 叶 变 换 组 的 个 数 。 
经 过 上 述 时 延 改 正 后 ， 再 进行 二 次 项 时 延 系 数 的 拟 合 。 设 未 知 的 二 次 项 时 延 系 数 为 
[Xo x1 xz]， 求 解 时 延 的 多 项 式 为 
a” = xo + xiti + x,t? (15) 
其 中 ，ti =ixT, 为 每 个 二 维 传 里 叶 变 换 组 开始 的 时 间 。 
对 上 式 右 端 求 导 ， 得 到 时 延 率 的 多 项 式 为 
bi = x1 + 2x2; (16) 
Wear sk (15) 和 式 (16) 建立 超 定 方程 组 ， 再 利用 QR 分 解 求解 ”。 


2 算法 验证 


RadioAstron 项 目 是 苏联 在 20 世 纪 80 年 代 中 期 提出 的 一 项 空间 其 长 基线 干涉 测量 
(SVLBI) 计划 。2011 年 ， 俄 罗斯 发 射 了 10m 口 径 的 RadioAstron 空 间 射电 望远镜 (Spektr-R, 
Ra) FÆ .Ra 与 地 面 多 个 大 型 射电 望远镜 构成 了 空间 VLBI 系 统 。 其 接收 频段 为 0. 3GHz、1. 6GHz、 
5. 06Hz 与 22GHz， 轨 道 高 度 近地点 600km， 远 地 点 330000km， 通 过 联合 地 面 望远镜 进行 干涉 测 
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为 验证 空间 VLBI 射 电源 条 纹 搜 索 算 法 ， 我 们 使 用 了 2014 年 空间 射电 望远镜 Ra 与 美国 阿 
雷 西 博 射电 望远镜 (Arecibo Radio Telescope, Ar) 联合 观测 射电 源 0823+033 的 数据 进行 
干涉 处 理 。 观 测 起 始 时 间 为 2014 年 4 月 11 日 23 时 07 分 00 秒 ,结束 时 间 为 2014 年 4 月 11 日 23 时 07 
分 59 秒 。 其 中 Ar 站 为 2 比特 量化 ，4 比 特 肩 出 的 数据 ， 衬 间 站 Ra 为 1 比特 量化 ，2 比 特 扇 出 的 数 
据 。 天 衬 频 率 4836MHz， 带 宽 16MHz， 积 分 时 间 为 8.192 x 10 习 秒 ， 相 关 处 理 后 的 数据 规格 为 
60 x 120 x 8192， 利 用 每 120 x 8192 个 点 的 数据 为 一 组 进行 二 维 傅 里 叶 变 换 。 由 式 〈10) 
ASK (11) 确定 时 延 和 时 延 率 方向 的 搜索 的 步 长 为 : 


N 一 4 
dstep = gg = 1.02400 x 107s 
M -8 
drstep = ZAT = 1.26210 x 10 s/s 
0 


以 先 验 时 延 模型 作为 搜索 的 中 心 点 ， 时 延 和 时 延 率 搜索 的 网 格 大 小 为 20 x 200， 其 中 时 
延 的 搜索 范围 为 : 6.33600 x 10-3s~8.28160 x 10-3s， 时 延 率 的 搜索 范围 为 : 1.62700 x 
10-6s/S~4.13858 x 10-6s/s。 

算法 验证 时 ， 首 先 采用 常规 多 重 网 格 条 纹 搜索 算法 ， 仪 使 用 线性 时 延 模型 ， 对 每 个 网 
格 点 相关 处 理 , 得 到 峰值 点 的 干涉 条 纹 图 。 再 使 用 本 文 研 究 的 射电 源 条 纹 搜索 算法 构建 二 次 
项 时 延 模型 ， 最 后 对 两 者 的 结果 比较 验证 。 

图 4(a) 为 直接 利用 线性 时 延 模型 进行 条 纹 搜索 的 相关 幅度 图 ， 其 峰值 点 对 应 的 网 格 位 
HY (15, 158) ， 相 关 幅 度 值 为 2.86900 x 1074, Bl4) 为 其 峰值 点 对 应 的 条 纹 相 位 。 
中 没有 明显 干涉 条 纹 。 
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(a) (b) 
图 4 使 用 条 纹 搜 索 和 线性 时 延 无 法 得 到 正确 条 纹 。 (a〉 相 关 幅 度 图 ;，〈b) 干涉 条 纹 图 


Fig. 4 Fringe search using linear delay can not produce the fringe. (a) Correlation amplitude; (b) interference fringe 


图 5 为 空间 VLBI 射 电源 条 纹 搜 索 算法 得 到 的 不 同 网 格 点 对 应 的 相关 幅度 值 
应 的 网 格 位置 为 11，101〉， 相 关 幅 度 值 为 1.96700 x 10-3， 经 过 拟 合 后 得 到 二 次 项 时 延 
系数 [7.36875 x 10-3s，2.88784 x 10-6s/s， 一 4.30446 x 10-11s/s?|。 使 用 该 二 次 项 模型 
做 相关 处 理 ， 可 得 到 图 6 (a) 的 干涉 条 纹 。 图 6 (b) 为 解 卷 绕 之 后 的 条 纹 结果 。 可 见 ， 本 文 提出 
的 算法 在 处 理 RadioAstron 数 据 时 要 优 于 多 重 网 格 条 纹 搜索 算法 。 
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Fig.6 Clear 
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图 5 二 次 项 时 延 补偿 的 相关 幅度 


Fig.5 Corre 
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(a) 
用 二 次 项 时 延 补偿 后 获得 清晰 的 干涉 条 纹 。 


该 算法 原理 样机 在 Linux 集 群 


unwrapping 


平台 ， 采 用 计算 节点 间 MPI+ 节 点 内 POSIX 多 线程 的 两 级 并 
行 实现 ， 计 算 节 点 的 CPU 型 号 为 英特尔 的 Xeon E5-2620， 主 频 为 26Hz， 每 个 节点 的 CPU 个 数 为 


Delay/s 


ChinaXiv 合 作 期 刊 


ation amplitude using quadratic delay 


As-Ra fringe phase - 


(b) 


(a) 解 卷 绕 之 前 的 干涉 条 纹 ;，(b) 解 卷 绕 之 后 的 干涉 条 纹 


interference fringes are obtained by quadratic time delay compensation. (a)Phase wrapping; (b)Phase 


12， 运 行 耗 时 如 表 1 所 示 。 后 续 通 过 优化 ， 运 行 速度 可 进一步 提升 。 


表 1 空间 VLBI 射电 源 条 纹 搜索 算法 运行 耗 时 


Tab. 1 Time consuming of the radio source fringe search algorithm for space VLBI 


Number of nodes*Number of threads 


Grid points (delay*delay rate) 


Time consuming 


10*10 


20*200 


About 45 minutes 
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在 将 来 探 月 四 期 任务 中 ， 地 月 空间 YLBI 试 验 系统 受 中 继 星 轨道 扰动 和 设备 时 延 影响 ， 

可 能 会 导致 射电 源 预 报时 延 模 型 和 实际 时 延 有 较 大 差异 。 本文 在 常规 多 重 网 格 条 纹 搜索 基础 
E, 根据 空间 VLBI 信 号 弱 的 特点 ， 提 出 了 一 种 适用 于 空间 VLBI 的 射电 源 条 纹 搜索 算法 ,该 算 
法 先 建 立 关 于 时 延 和 时 延 率 的 二 维 搜索 网 格 , 使 用 线性 时 延 相 关 处 理 , 然后 由 二 维 傅 里 叶 变 
换 搜 索 残 余 时 延 和 时 延 率 , 最 后 重 构 二 次 项 时 延 模型 。 本 文 使 用 RadioAstron 数 据 验 证 算法 ， 
结果 表明 在 使 用 常规 线性 时 延 模型 搜索 无 法 实现 空地 基线 条 纹 的 情况 下 , 新 算法 可 以 得 到 清 
晰 的 空间 VLBI 射 电源 干涉 条 纹 。 经 进一步 测试 和 完善 ， 该 算法 可 应 用 于 我 国 后 续 空间 VLBI 
项 目 中 。 
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Abstract: The fourth phase of China's lunar exploration project will build the first lunar 
terrestrial Space VLBI test system composed of 4.1m antenna on Chang'e 7 relay satellite and 
ground-based radio telescope. Different from the ground-based VLBI system, due to the influence 
of relay satellite orbit disturbance and equipment delay, the predicted delay model of radio source 
may be quite different from the actual delay, which can not guide the relevant processors as the 
core equipment of VLBI data processing to work normally. In this paper, a radio source fringe 
search algorithm for Space VLBI is proposed. Aiming at the characteristics of weak signal 
receiving ability of space small aperture antenna and the need to use quadratic delay model to 
ensure long-time integration accuracy, firstly, a two-dimensional search grid about delay and 
delay rate is established, Using the prior delay and delay rate information of each grid point, the 
linear delay is constructed to process the original data in different time periods, and then the 
residual delay and delay rate are searched by two-dimensional Fourier transform to determine the 
grid point corresponding to the maximum correlation amplitude value. Finally, the quadratic delay 
model is reconstructed by using the delay corresponding to the grid point. The algorithm can be 
applied to the radio source fringe search of the ground air baseline of the RadioAstron project. 


Key words:Space VLBI, Fringe search; Delay model; Correlation processing; Radio source 


